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^ , Abstract 

■ We consider a dynamical system for computing Nash bargaining solutions on graphs 

and focus on its rate of convergence. More precisely, we analyze the edge-balanced dy- 



^ ■ are quadratic in the number of matched edges. 
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namical system by Azar et al and fully specify its convergence for an important class of 
elementary graph structures that arise in Kleinberg and Tardos' procedure for comput- 
ing a Nash bargaining solution on general graphs. We show that all these dynamical 
systems are either linear or eventually become linear and that their convergence times 



1 Introduction 

Bargaining and, in particular, the concept of Nash bargaining on general graphs has been the 
focus of much recent research in economics, sociology and computer science [21 [51 [IHl El [12] • 
^ ■ In a bargaining system, players aim at making pairwise agreements to share a fixed wealth 
O . specific to each pair of players. Bargaining solutions provide predictions on how the wealth 

will be shared and how this sharing would depend on players' positions in a network describing 
some notion of relationships among players. 

The concept of Nash bargaining solution was introduced by Nash |T3] for two players, each 
having an exogenous, alternative profit at its disposal were they to disagree on how to share 
the wealth. Recent research has focused on the concept of Nash bargaining with multiple 
players where each player has alternative profits determined by trading opportunities with 
neighbors in a graph. In the computer science literature, Kleinberg and Tardos [12j were the 
first to establish various properties of Nash bargaining outcomes on general graphs. They 
also propose a polynomial-time algorithm for computing them, provided one exists. Follow 
up work aimed at introducing some local dynamics that are natural (so they, hopefully, have 
some connections with reality) and studied their convergence properties. In particular, Azar 
et al [2] considered the so called edge-balanced dynamics and established various properties 
about fixed points and convergence but left open the characterization of the convergence rate. 
In a tandem of papers jHllin], Kanoria et al considered an alternative, natural dynamics, and 
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established polynomial convergence time bounds under some generic assumptions. An open 
research question has been to gain a better understanding of convergence properties and 
obtain tight bounds on the convergence time for these types of systems. 

In this paper, we consider the edge-balanced dynamics of Azar et al ||2j over elementary 
graphs that arise in the decomposition procedure of Kleinberg and Tardos |12] which include 
a path, a cycle, a blossom and a bicycle (see Figures [H [21 E] and IHfor examples). It turns 
out that, for all these network structures, the dynamics is either linear or eventually becomes 
linear. Specifically, we show that the dynamics is linear for a path and a cycle and is 
eventually linear for a blossom and a bicycle (and characterize the time when this takes 
place). This allows us to fully characterize the rate of convergence by deploying well known 
spectral methods for linear systems. As a result, for all these elementary structures, we find 
that the convergence time is quadratic in the number of matched edges. 

1.1 Outline of the Paper 

In Section [2] we introduce system assumptions and overview relevant concepts, including the 
concept of Nash bargaining outcomes, local dynamics, and the KT procedure. Section [3] 
provides the characterization of the edge-balanced dynamics and convergence times for each 
of the elementary graphs of the KT decomposition. Finally, we discuss related work and 
conclude in Sections H] and |5] respectively. 

2 System and Assumptions 

2.1 Nash Bargaining Outcomes on Graphs 

We consider a graph G = (V, E) where V is the set of nodes and E is the set of edges. Each 
node corresponds to a distinct player that participates in the trading game defined as follows. 
Each edge {i,j) & E is associated with a weight Wij > representing the amount that can 
be shared between players i and j should these two players decide to trade with each other. 
The trading game is one-exchange meaning that each player attempts to make a pairwise 
agreement with at most one other player, which corresponds to a matching M G E in the 
graph where (z, j) G M if and only if players i and j reached an agreement. We denote with 
Xi the profit of player i where Xj > and let x = (xj)igy denote the vector of players' profits. 

A balanced outcome or a Nash bargaining solution is a pair (M, x) where M is a matching 
in G and x is a vector of players' profits. Such an outcome satisfies the following properties: 

• Stability: for every edge {i,j) G E, 

Xi -\- Xj ^ 

• Balance: for every {i,j) G M, it holds that 

Xi - max {wi^k - Xk)+ = Xj - max (w,- ^ - Xk)+ 

where, hereinafter, Vj denotes the set of neighbors of a node i and (•)+ := max(0, ■). 
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The stability property means that there exists no player that can improve her profit 
by unilaterally deciding to trade with an alternative trading partner. The balance property 
originates from the Nash bargaining problem [14j where two players 1 and 2 aim at a pairwise 
agreement to share a profit w having outside profit options ri and r2 in case of disagreement. 
The Nash bargaining solution is then for players 1 and 2 to share the surplus w — ri — r2 
equally, if positive, i.e. receive profits Pi = ri + \{w — ri — r2)+ and P2 = ^2 + \{w — ri — r2)+, 
respectively. This allocation is balanced in the sense that Pi — ri = p2 ~ "^2) which is exactly 
the above asserted balance property where the outside profit options are determined by the 
values that players may extract through trading agreements with their neighbors. 

2.2 KT Procedure 

Nash's bargaining solutions on graphs are intimately related to maximum-weight matchings. 
In [12j it was found that the matching M of a stable outcome x is a maximum-weight 
matching. Furthermore, whenever a stable outcome exists, a balanced outcome exists as 
well [12j . The outcome vector x can be seen as a feasible solution of a dual to the fractional 
relaxation of a maximum- weight matching (primal): 

maximize I](ij)ei7; WijXi^j 

over Xij > 0, [i,]) G E 
subject to Ei:(ij)ei?^M ^ 1 

where the dual problem is the following linear problem with two variables per inequality: 

minimize Xliev 

over Xj > 0, i & V 
subject to Xi + Xj > Wij, G E. 

In [T2j , it was established that a balanced outcome (M, x) can be found in polynomial 
time by first finding a maximum- weight matching M and then solving the above dual problem 
to find a balanced vector x. The dual problem can be solved by an iterative procedure where 
each iteration maximizes the smallest slack as described in the following. 

We denote by Si the slack of node i defined by 

Si = Xi- max (wj i -xi)+ 
{i,i)eE\M 

while the slack of edge denoted by Sij is defined by 



Indeed, for a stable outcome x, Si,j > 0, for every € E. It is not difficult to observe 

that node and edge slacks satisfy 

Si = mm{xi, min Si^i) . 

{i,l)&E\M 

The KT procedure for finding a balanced outcome proceeds by successively fixing the 
values Xi for some nodes in V . This is allowed by the following key property |12j: if there 
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exists a set y4 C and a > such that Si < a for every i & A and Si > a for i ^ V \ A 
and a vector x such that values Xi are balanced in A, then there exists a vector x' such that 
= Xi for every i & A that is a balanced outcome for G. 

The KT algorithm is sketched as follows. Let a > be a variable and let A be a set of 
nodes for which values Xi have been already assigned. The set A is constructed such that 
no matched edge crosses the cut {A, V \ A), i.e. for every node i & A there exists no node 
j & V \A such that G M. Initially, cr = and the set A contains all the unmatched 
nodes. The algorithm then proceeds inductively with respect to the number of nodes with 
unassigned values as given by Given a and A the inductive step amounts to assigning 

values to nodes in V\A that maximize the minimum slack a' > a which amounts to solving 
the following linear program 

maximize a' 

subject to x[> a', i E V \A 

x'i + x'j = Wij, e M (1) 

x'i + x'j > w,j + a', (z, j) G E \ {M U E{A)) 

where £"(^4) corresponds to the set of edges of the graph G linking nodes in A. 

For a fixed o"', this is a linear inequalities' problem with at most two variables per inequal- 
ity, for which polynomial algorithms exist. In particular, by results of Aspvall and Shilach [1], 
for a given o"', the system of inequalities is infeasible if there exists an infeasible simple loop 
in the graph construction described in [Ij. A path is said to be a loop if the initial and final 
nodes are identical and is said to be simple if all intermediate nodes of this path are distinct. 
Furthermore, if a feasible solution exists than it can be constructed by finding the most con- 
straining feasible simple loop. For the above system of inequalities, any such feasible simple 
loop is either a path, a cycle, a blossom or a bicycle. We refer to these as KT elementary 
graphs and define them in the following: 

• Path. A path consists of alternating matchings with each of its end nodes anchored at 
either a node z G A or at a matched edge («, j) G M such that Sj = Xj^ 

• Cycle. A cycle consists of an even number of nodes connected by a path of alternating 
matchings. 

• Blossom. A blossom is a concatenation of a cycle and a path (we refer to it as a stem) 
as follows. The cycle consists of an odd number of nodes that are connected by a cycle 
of alternating matchings started from a node (we call gateway) with an unmatched 
edge. The stem is a path of alternating matchings such that one end node is matched 
to the gateway node and the other end node is anchored as for a path. 

• Bicycle. A bicycle is a concatenation of two blossoms by connecting the end nodes of 
their respective stems such that the connected stems form alternating matchings. 

The above described step is repeated until all the nodes are assigned values, i.e. until 
V \ A = 0. Hence, the total number of such steps k is at most the number of nodes \V\. At 

^Recall that if for a matched edge E, i ^ V \ A, then also j ^ V \A, and vice versa. 
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each step /, a KT elementary structure Ci and maximum slack ai are identified such that the 
(T;'s form an non-decreasing sequence, = (Jq < cxi < ■ ■ ■ < ak- 

2.3 Convergence 

We introduce a few elementary concepts about stability of dynamical systems in a somewhat 
informal manner and then define the notion of convergence time considered in this paper. 
We say that a dynamical system, according to which x{t) evolves over t > 0, is asymptotically 
stable, if for every initial value x{0), there exists a point if* such that 

lim \ \x{t) — x*\ \ = 0. 

The system is said to be globally asymptotically stable if x* is unique, i.e. does not depend 
on the initial value x(0). 

In particular, for a linear system 

x{t+l)= Ax{t) + hit) 

where A is a given matrix and h{t) is a vector that may depend on t, we have that the system 
is globally asymptotically stable if the spectral radius of the matrix A is smaller than 1 (i.e. 
all eigenvalues are of modulo strictly smaller than 1). The concepts of asymptotic stability 
and global asymptotic stability are standard, see [Hj for more details. 

We say that the convergence to a point x* is exponentially bounded if there exist C > 
and -R > such that for every initial value x(0), we have 

\\x{t) - < Ce'^\ for every t > 0, 

where we refer to R as the rate of convergence and call T = 1/R the convergence time. 
Moreover, If x{t) evolves according to the aforementioned linear system then the rate of 
convergence is given by (i) R = log(l/p(A)) where p(A) is the spectral radius of matrix A 
if the system is globally asymptotically stable, and (ii) R = log(l/A2(A)) where A2(A) is 
the modulus of the largest eigenvalue of matrix A that is smaller than 1, if the system is 
asymptotically stable. 



3 Edge-Balanced Dynamics for KT Elementary Graphs 



The edge-balanced dynamics was first considered by Rochford |15] and Cook and Yamag- 
ishi [6j, this dynamical process assumes that players already agreed on a matching M and 
are negotiating the value of the outcome x. Hence, each matched player i is assigned a trading 
partner, which we denote with pi. A version of this dynamics in discrete-time can be repre- 
sented as follows. For a fixed < a < 1 and given an initial value x{0), for z = 1, 2, . . . , 
and t = 0, 1, . . ., we have that 



Xiit + 1) = Xi{t) + a 



Viit) + -{wi^p^ - yi{t) - ypXt)) 



"''Pi 



Xiit) 



(2) 
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where yi{t) is the best alternate value that a matched player / may get at time t by trading 
with her other neighbors, i.e. 



yiit) 




max {wi^k - Xk(t)) + 
',k)eE\M 



and we use the notation [-j^ = min(max(-, a),b), for a < b. 

It is not difficult to observe that if players i and j are matched, then Xi{t) + Xj{t) = Wij 
is time invariant, i.e. if the latter holds for a time t, then it still holds for time t + 1. 
Note that the dynamics is not necessarily consistent with Nash bargaining solution for every 
time t as for a matched pair the edge-surplus Wij — yi(t) — Ujit) is allowed to be 

negative; the only requirement is that the allocation yi{t) + ^{wij — yi(t) —yj(t)) is in [0, Wij]. 
However, the edge surpluses are guaranteed to be positive for t large enough [2]. Besides the 
dynamics described in ([2]) is guaranteed to converge to a fixed point that corresponds to a 
Nash bargaining solution, see [2, Theorems 1,2]. 

In this section, we will observe that for all the elementary graphs of the KT decomposition, 
the values held by the nodes eventually evolve according to a linear discrete-time dynamical 
system, i.e., for a given matrix A and a vector b{t), x{t) evolves according to 



We will find that for a path and a cycle the dynamics is linear for every time t > while 
for a blossom and a bicycle there exists a finite time Tq > such that the dynamics is linear 
for every t > Tq. The asymptotic behavior is determined by spectral properties of matrix A. 
Note that it suffices to consider the spectrum of matrix A for a = 1. This is because, for 
every given 0<a<l, A' = l— a + a A is an eigenvalue and v is an eigenvector of the matrix 
A, where A is an eigenvalue and v is an eigenvector of the matrix A under a = 1. We will see 
that for every KT elementary graph, the eigenvalues of matrix A, under a = 1, are located 
in the interval [—1,1] and will show that —1 can be an eigenvalue only for a cycle with an 
even number of matched edges or a bicycle with an even number of matched edges in each 
of its loops. In the latter two cases, for a = 1, there is no convergence to a limit point as the 
asymptotic behavior is periodic because of the eigenvalue —1. This is ruled out by choosing 
the smoothing parameter < a < 1, making all the eigenvalues strictly larger than —1, and 
thus ensuring convergence to a limit point. 

It is also worth noting that if the system is globally asymptotically stable and if when, 
the term b{t) in Q is equal to 6 a constant, that does not depend on t, then the fixed point 
to which the system converges is given by {I — A)~^b where {I — A)^^ = X]a:>o would 
like to mention that the analysis of the fixed points of these dynamical systems is relevant 
to understanding the bargaining power of nodes depending on their position in the network. 

Finally, we note that for our results in this section, we assume uniform edge weights and 
under this assumption, without loss of generality, we let We = 1, for every e E E. 



We consider a path with n matched edges with boundary values x~^{t) and x (t) as illustrated 
in Figure [TJ In this case, the evolution of the node values x{t) boils down to a discrete-time 



x{t+ 1) = Ax{t) + b{t). 



(3) 



3.1 Path 
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linear dynamical system 



where A is the n x n symmetric tridiagonal matrix 
/ 1/2 ••• \ 

1/2 ••. ••. •-. ; 

A= ■•• ■•• ■•• 



V 



••• •-. 1/2 

1/2 y 



(4) 



and b{t) = C-^,0,...,0,q^f 



n-2 



Figure 1: A path with boundary conditions. Note that the values at the ends of each matched 
edge are given by Xi and 1 — Xi, the latter being omitted in the figure. 



The eigenvalues of matrix A are = cos (;^) , k = 1,2, . . . ,n, with the corresponding 
orthonormal eigenvectors 

Note that every eigenvalue is of modulo smaller than 1. This implies asymptotic stability for 
every < a < 1. From the above spectrum, we have the following characterization of the 
convergence time: 

Theorem 1 For a path of n matched edges and every < a < 1, the convergence time is 

T = ^n''-[l + Oil/n% 



Proof. An eigenvalue A and associated eigenvector v of matrix A satisfy 

Xvi = ^V2 

Xvi = -Vi_i + -^Vi+i, l<i<n 

Using A = cos(0) and Vi = sm{(f)i), for > in the above equations, along with some 
elementary trigonometric calculus, it readily follows that = for k = 1,2, . . . ,n. 

Since — 1 < A^ < 1 for every k and Ai > has the largest modulo, the convergence time 
is given by T = log(l — a + aAi). Noting that Ai = 1 — ^ + 0(l/n^), the asserted result 
follows. ■ 
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From this theorem, we observe that the convergence time is quadratic in the number of 
matched edges. Moreover, if x+ = x_ = then it is not difficult to see that the dynamics 
convergence to the fixed point given by 

(/-A)-i(l/2,0,...,0f. 



3.2 Cycle 

For an ahernating cycle between n nodes (n even), the dynamics of node values x{t) boils 
down to a linear dynamical system ([2]) where A is the following circulant matrix, for n = 2, 
^,01 



1 



and otherwise 



A 



/ 

1/2 






V 1/2 



1/2 


1/2 





1/2 







1/2 \ 





1/2 




1/2 


1/2 





1/2 






and where vector h = 0. Note that in this case x{t) = A*x(0), for t > 0. 



(5) 




Figure 2: A cycle. 

By using similar arguments as for a path, it is not difficult to establish that the eigenvalues 
of matrix A are = cos ^^^^^^-^j , k = 1,2, ... ,n, with the corresponding orthonormal 
eigenvectors 

^(l,l,...,l,lf , if A;=l 
v, = <^ ^_(l,-l,...,l,-lf , ifk = l+n/2 

(1, cos (0fc) , . . . , COS {(f)k{n - 1))) , o.w. 



where for ease of notation, 0^ = ^^(^ ^)(" ^) _ 



Using the spectral decomposition of the symmetric matrix A (see [TT] for details), we 
have 

n 

xit) = j2>^ivkvim- (6) 

k=l 



We distinguish two cases: 

• Case 1: n is even. In this case, = —1, for k = 1 + n/2, and > —1, otherwise. 
From ([6]), we have 

X{t) = [viV^ + (-l)Vi+„/2tTf+„/2) x(0) + o(l). 

Therefore, the asymptotic behavior is periodic. 

• Case 2: n is odd. In this case, —1 < A^ < 1, for every k, and thus we have asymptotic 
convergence to the limit point, \imt^ooXi(t) = ^ Yl]=i^ji^)y every i. 

In view of the above observations, we note that for even n, we need to assume that a is 
strictly smaller than 1 in order to rule out asymptotically periodic behavior, while for odd 
n, we can allow for a = 1. Moreover, it is worth noting that the limit point, in both cases, 
depends on the initial condition and is given by ^ Yl^=i • • • 5 1)"^- 

The following result shows that in like manner as for a path, the convergence time is 
quadratic in the number of matched edges, but note that it is four times smaller, asymptot- 
ically for large n. 

Theorem 2 For cycle graph with n matched edges and for a G (0, 1) the convergence time 
is 

T=^^n'-[l + 0{l/n% 

3.3 Blossom 

A blossom is a concatenation of a cycle and a path (we refer to it as a stem); see Figure [3] 
for an example. We consider a blossom with n matched edges in the stem and m matched 
edges in the loop. We refer to the node that connects the stem and the loop as a gateway 
node. The matched edges of the stem are enumerated as 1,2, ... ,n along the stem towards 
the gateway note. We let Xi denote the value of the end node of an edge i of the stem that 
is connected to a node towards the open end of the stem. Similarly, we enumerate matched 
edges of the loop as 1,2, ... ,m and let yi denote the value of the node that appears first on 
a matched edge i as we go along the loop in the clockwise direction. 




Figure 3: A blossom. 
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It can be observed that node values x{t) and y{t) evolve according to the following non- 
linear dynamical system: 



Xi{t+l) 



X2{t) 

2 



Xi[t + 1) = ^ , 1 <i <n 



Xn{t + 1) 



1 + Xn-lit) - max[l - yi{t),y^{t)] 

^ (7) 

Xn{t)+y2{t) ^ ^ 

2 



yi{t + l) = , l<t<m 



Vmit+l] 



1 + yrn-lit) - Xnit) 



2 

Note that the system is non-linear only because of the maximum operator that acts in 
the update for the node matched to the gateway node, which connects the stem and the 
loop of the blossom. The maximum operator is over the values of the nodes that are in the 
loop and are matched to the neighbors of the gateway node, 1 — yi{t) and ym{t)- It turns 
out that, eventually, one of these two values is larger or equal to the other and, hence, the 
system dynamics becomes linear. This is shown in the following lemma. Note that the sum 
Viit) + limit)-! if smaller than or equal to 1, indicates max(l — ?/i(t), ymif)) = 1 ~ Viif)^ and 
otherwise, max(l - yi{t),ym{t)) = ym{t)- 

Theorem 3 For a blossom with n matched edges in the stem and m matched edges in the 
loop, for every initial value (x(0),y(0)), the sum of node values yi{t) + ym{t) satisfies: 

1. yi{t) + ym{t), for t > 0, is autonomous of x{t), t > 0. 

2. \imt^ooyi{t) +ym(t) = 1. 

2 

3. The asymptotic rate of convergence is . 

4- There exists a time Tq > such that either yi{t) + ym{t) < 1 or yi{t) + ym{t) > 1 for 
every t > Tq. 

5. Tn = Oim'^). 



As an aside remark, note that the value sum Xll^i f*^^ \oo^ nodes evolves au- 
tonomously from t > 0. To see this, from ([7]) note 

Y^y^{t + i) = Y, y^it) - 2 (2/1 W + ymit)) + 2- 

i=l i=l 

and that the claim follows from Theorem [3] item 1, saying that yi{t) + ym{t) evolves au- 
tonomously of a;(t), t > 0. 

The theorem derives from an explicit characterization of yi{t) + ym{t), for every t > 0, 
which we present in the following: 
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Lemma 1 Given initial value y{0), for every t > 0, 

r — 

I 2 



viit) + vut) = 1 - Yl f2i-im))A^i (8) 

m + 1 ^ — ^ 

j=i 



Proof. The part of the system y{t) evolves as the following non-autonomous linear system 
y{t + 1) = Ay{t) + b{t) where A is a tridiagonal matrix that corresponds to a path of m 
matched edges and b{t) = (a;„(t)/2, 0, . . . , 0, (1 - x„(t))/2)^. 

m-2 

Since A is a symmetric matrix, we can use the spectral decomposition A = YlT=i ^k^kV^ 
where Ai, A2, . . . , Am are the eigenvalues and Vi,V2, . . . ,Vm are the orthonormal eigenvectors 
of matrix A, which we identified in Section 13.11 

Using the spectral decomposition, we note 

m t—l 
k=l s=0 



where v^^i = y sin {:^^^i) is the i-th coordinate of the eigenvector Vk- Summing up yi{t) 
and ymit)) we obtain 

m / t-i 

yi{t) + ym{t) = J2^Vk,l + Vk,m) Kv^m + J2 >^V~^^lKs) 
k=l \ s=0 

k odd V .=0 

where the last inequality is because of the fact Vk,m = Vk,i for k odd and Vk,m = —Vk,i for k 
even. Furthermore, 

-rr/ X Xnis) 1 - Xn{s) Vk,l 

v^b{s) = Vk,i— \-Vk,m — lor k odd. 

Therefore, 

yi{t)+yut) = Y (^>k,iv^m + Y^i~'~'<i 

k odd V *=o 



Y U>kAm) + yY<') 

fc odd ^ ^ ^ 

k odd ^ 



2 



l-A2Cy(0) + (l-Ai,)^(l + A, 

tY(' + ^^)-^,Y (l + A.-2y^v^^ty(0)]A*. 
k odd k odd \ / 
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It remains only to show that 

^ E (i + ^^) = i 

m + 1 ^-^ 

k odd 

which follows readily by elementary trigonometric calculations. ■ 

We now give the proof of Thereom |3l 

Proof of Thereom [31 The statements of the theorem derive from Lemma [T] as follows. 
Item 1 clearly holds as the function ([H]) depends only on the initial value y{0). Item 2 
follows from (|8]) because all the eigenvalues are real and with modulo strictly smaller 
than 1. Item 3 holds from the fact that the largest modulo eigenvalue of matrix A is Ai = 
cos {^) = 1 - i + 0{l/m^) and hence R = log(l/Ai) = ^ + 0{l/m^). Item 4 holds as 
the sum in ([8]) is asymptotically dominated by the largest modulo eigenvalue X2i-i such that 
v^i-iVi^) 7^ 0, i.e. the mode associated to the eigenvalue X2i-i is excited. Let us consider 
the case where such an eigenvalue is Ai and m is even; the other cases follow by similar 
arguments. From Lemma [T|, we have yi{t) + ymit) = 

([—1 
i=2 



m + 




and, thus, since |A2i-i/Ai| < 1, for every 1 <i < [m/2], for t large, 

yi{t) + ym{t) = 1 - ^rr^i [/i(y(0)) + 0(1)] 

m + 1 

Finally, item 5 holds as, for m large enough. 



|A2i-i| A3 An' 4 

7 := max — — < -;— = 1 7 + (J[l/m ). 

Ai Ai 



For an arbitrary e > 0, we have |A2t-i/Ai|* < e, for every i > 1, provided that time t is such 
that 



47^2 



t > = -izT^rn'll + 0(1)]. 

log 



. 7 

Hence, Tq = 0{rn?). 



From Theorem [3] item 4, we have that the dynamics for a blossom is eventually according 
to the following linear system 

y{t + i) )-^\yit))^^ 

where matrix A and vector h assume one of the following two choices: 
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A 



(9) 



Case 1: (1 -yi(t) > y^{t)) 

T„ P 
Q T„ 

with T„ and tridiagonal matrices of paths of n and m matched edges, respectively, 
and 

0\ /O 



/ 




v 




\ 










Q 



v 










\ \ 






-i 



and h= (0^_^^, 1/2)^. 

n+m—l 

Moreover, if the initial condition is such that 1 — ?/i(0) > ?/m(0), e.g. yi{t) = for all i, 
then the fixed point is given by (/ — A)^^b. 

Case 2: (1 — < ym{t)) same as under Case 1 but 



/ 





V 



\ 






-i 



and b = (0, . . . , 0, 1/2, 0, . . . , 0, 1/2)^. Finally, it the initial condition is such that 



n-l 



m—l 



1 — yi{0) < ym{0), e.g. yi{t) = 1 for all i, then the fixed point is given by (J — A) ^b. 

In the following we only consider Case 1 as the spectrum of matrix A under Case 2 is ex- 
actly the same. We note that the eigenvalues of the matrix A are (Ai, A2, . . . , Xn+[m/2i , /^i, f^2, • • • , ] 
where 



A. 



cos 



COS 



27Tk 



2n + m + 1 
7i{2k - 1 



m + 1 



, k = l, 
k = l,.. 



. , n + [m/2\ 
\m/2] 



with a proof provided in Section 16.11 

It is noteworthy that all the eigenvalues have modulo strictly smaller and 1, and thus, the 
system is globally asymptotically stable. We now characterize the convergence time from an 
instance at which the system became linear. 

Theorem 4 For a blossom with n matched edges in the stem and m matched edges in the 
loop, for every < a < 1, the convergence time T satisfies: if m is even, then 



:{2n + my ■ [1 + 0(1)] 
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otherwise, for m odd, 



T = — - max ^m^, ^{2n + mf^ ■ [l + o(l)]. 

Observations. The result implies that the convergence time is 0((n+m)^), i.e. quadratic 
in the number of matched edges. There is a significant difference with regard to whether the 
number of matched edges in the loop, m, is even or odd. The convergence is slower for m even. 
Specifically, if the length of the stem is at least twice the length of the loop, the convergence 
time is larger for a factor 4. For a fixed n, the convergence time is asymptotically -^^"^ as for 
a path of length m which is intuitive. Likewise, if m is fixed and odd, the convergence time 
is asymptotically -^n^ as for a path of length n and thus also in conformance to intuition. 

Proof. For the eigenvalues Ai, A2, . . . , )^n+[m/2i , it is readily checked that 

, , /(I + 1 r^AA )vr' 

max Afc — -An+lm/2\ — COS — - — ■ — 

k \ 2n + m + 1 

where 1^ stands for the indicator that m is odd, while, on the other hand, 

TT 



max = /ii = cos 
k \m + 1^ 

Therefore, the spectral radius of matrix A, p(A) = max(maxfc |Afc|,maxfc \nk\) is given by 

f cos (2;^:^), meven 
p(A) = <^ cos i ^n+^+i ) > m odd , m < 2n - 1 
cos {-;^^) , odd ,m > 2n — 1. 

The asserted asymptotic follows from the last identities. 



3.4 Bicycle 



>y2 



\3 Zl '° 

Figure 4: A bicycle. 

A bicycle graph consists of two loops that are connected by a path. Without loss of 
generality, we refer to one of the loops as loop 1 and to other as loop 2 and refer to the 
path as a cross-bar; see Figure H] for an illustration. Notice that a bicycle graph corresponds 
to a concatenation of two blossoms by connecting the end nodes of their respective stems 
so that a cross-bar is formed of alternating matchings. We let / and m be the number of 
matched edges in loop 1 and loop 2, respectively, and let n be the number of matched edges 
of the cross-bar. The values of the end nodes of the matched edges are denoted by z{t) = 
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{zi{t),Z2{t),...,zi{t))'^, x{t) = {xi{t),X2{t),. . .,Xn{t))'^ and y = {yi{t),y2{t), . . . ,ym{t)V for 
loop 1, cross-bar, and loop 2, respectively. See Figure H] for positions of the corresponding 
nodes. 

We note that for a bicycle the system evolves according to the following non-linear system: 

.Kt + i) = ^^-i^^^±^ (10) 



+ 1) 



2 

_ X2{t) + max[l - zi{t), zi{t)] 
~ 2 

plus other updates as for blossom (j7]). 



In this case, the non-linearity originates because of two gateway nodes that connect the cross- 
bar with loops, each such gateway node having two alternative profit options with nodes in 
the loops. Similarly as for a blossom we have that eventually the dynamics becomes linear 
as stated in the following: 

Proposition 1 For a bicycle with I and m matched edges in loops and n matched edges in 
the cross-bar, there exists a time Tq > such that for every t > To, {z{t),x{t),y{t)) evolves 
according to a linear system. Furthermore, Tq = 0(max(/, m)^). 

This observation follows from Theorem |3] applied to each loop of the bicycle. This can be 
done because both yi{t) +ym{t) and Zi{t) + Zi{t) evolve autonomously as given by Lemma [1] 
for yi{t) + ymit) and analogously for zi(t) + zi{t). 

We have shown that the dynamics for a bicycle is eventually according to a linear system, 
which is specified as follows: 



fill 





where 

A = 

with the given matrix blocks defined by 





Q' 





P' 




p 





Q 


T 



Ti Q' ^ / T„ P 



are the matrices that correspond to two blossoms formed by loop 1 and cross-bar, and cross- 
bar and loop 2, respectively. 

The pair (A, b) admits four possible values, corresponding to all possible combinations of 
two cases for each of the loops (Case 1 and Case 2 in Section |33|) : 
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1. Both P' and P as in Case 1 
6 = (l/2,0^_^,l/2f, 

l+n+m—2 

2. P' as in Case 1, P as in Case 2, 

h= (l/2,0^_^,l/2,0^_^,l/2f , 

l-l m+n-2 

3. P' as in Case 2, P as in Case 1, 

h= (l/2,0^_^,-l/2,0^_^,l/2f , 

l+n—2 m—1 

4. Both P' and P as in Case 2, 

b= (l/2,0^_^,l/2,0^_^,-l/2,0^_^,l/2f . 

l-l ra-2 m-1 

In the following, we will only consider the case under item 1 as the same end results hold 
for other cases. The eigenvalues of the matrix A are given by 



/7r(2^-l^\ A; = l,...,[m/21, 
\ m + 1 / 

'=°K^^rflT^)-'= = ^--"+L'/2J + L.n/2J 

which we establish in Section 16.21 

Remark that in any case all the eigenvalues are strictly smaller than 1. On the other hand, 
if both / and m are even, then —1 is an eigenvalue with eigenvector (1, —1, 1, —1, . . . , 1, —1)"^, 
and otherwise, all the eigenvalues are strictly larger than —1. Therefore, if both / and m are 
even, then the asymptotic behavior of system (1111) is periodic, while otherwise, it is globally 
asymptotically stable. 

As a byproduct, similarly to Theorem |U we can establish that from an instance at which 
the system became linear, the convergence time scales as follows. 

Theorem 5 For a bicycle with n matched edges in the stem and m and I matched edges in 
the loops, the convergence time T satisfies the following. If m or I is even, then for every 
< a <1, 

2 



T = ^(2n + m + /)2- [l + o(l)] 
otherwise, if m and n are odd, then for every < a < 1, 



T = ^max ( m^/^-(2n + m + /)M ■ [l + o(l)]. 
a-K^ \ 4 / 

Therefore, the convergence time is 0((Z + n + m)^), i.e. quadratic in the number of 
matched edges. 
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4 Related Work 



The concept of balanced outcomes was introduced by Nash in |14j for the case of two players 
with exogenous profit options. This concept follows from a set of axioms and different axioms 
were subsequently considered; e.g. see [T^ . 

Kleinberg and Tardos |12] considered the concept of Nash bargaining solutions on graphs 
where profit options available to a player are not exogenously given but determined by her 
position in the graph. They established relations between stable and balanced outcomes and 
devised a polynomial time algorithm for computing balanced outcomes. Their work left open 
the question on existence and properties of local dynamics. 

A local dynamics for Nash bargaining on graphs was recently considered by Azar et 
al [2]. This paper assumed a fixed matching of nodes and considered a local, so called edge- 
balanced dynamics, for outcome vector x. They established that fixed points of this dynamics 
are balanced outcomes and that the convergence to the fixed point occurs in an exponential 
number of rounds. A concurrent and independent work by Cells, Devanur and Peres [4J 
considered the rate of convergence of edge-balanced dynamics. Their approach is different 
from ours in using a reduction of edge-balanced dynamics to a random-turn game for a class 
of graphs with uniform edge weights. For this class of graphs, their convergence time is 
quadratic in the maximum path length of an auxiliary graph derived from the input graph 
and given matching. This class of graphs does not accommodate cycles, but accommodates 
paths, blossoms and bicycles and for these cases the bound is quadratic in the number of 
matched edges. Another difference with our work is that for each of our elementary subgraphs 
we provide explicit characterization of dynamics and tight asymptotic estimate of convergence 
rates (exact constant factors). 

The assumption that matching is fixed was removed by Kanoria et al by introducing a 
natural dynamics studied in [H] and (TU]. It was shown in pj that provided that there exists a 
unique Nash bargaining solution and the graph satisfies the positive gap condition, the natural 
dynamics converges to this Nash bargaining solution in a polynomial time. Specifically, they 
showed that there exists a constant C > such that the convergence time is upper bounded 
by C[W/a + \og{a/e)]n^^^ , where W is an upper bound on the maximum edge weight, cr > 
is the gap and e,6 > 0. In [7j, we showed, using techniques similar to the ones introduced 
in this paper, that the bound in the number of nodes can be improved to 0{n^~^^). Using a 
different approach, in [lOj, the authors established that if the maximum weight matching is 
unique, then there exists T = 0{n^/g^) such that for every initial value the natural dynamics 
induces the maximum-weight matching, for every t > T; where g is the difference between 
the weight of the maximum- weight matching and the weight of the second best matching (we 
refer to as matching weight gap). 

Finally, another related work is that on maximum-weighted matchings on graph. There 
indeed is a close connection between stable outcomes and maximum weight matchings re- 
flected by the similarity of the distributed algorithms considered for solving both problems. 
In particular, Bayati et al [3] considered an auction-like algorithm, which is similar in spirit 
to the natural dynamics for solving the balanced allocation problem, and showed that for 
complete bipartite graphs with a unique maximum-weight matching, the convergence time is 
OiWn/ g) where W is the maximum edge weight, g is the matching weight gap and n is the 
number of nodes. 
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5 Conclusion 



In this paper we showed that some known Nash bargaining dynamics on graphs can (even- 
tually) be characterized by linear dynamical systems and this enabled us to derive tight 
characterizations of their convergence rates. Note that if the dynamics, as restricted to each 
of the KT elementary graphs that arise in the KT decomposition, were decoupled then the 
previous analysis will yield O(n^) convergence time, since there are at most n such structures 
each taking O(n^) time to converge. 

An interesting direction for future work is to investigate the extent by which the dynamics 
on the different KT substructures are coupled under assumptions such as the positive gap 
condition of the KT decomposition or the matching weight gap. Another interesting direction 
is to analyze the bargaining power of nodes based on their network position. 
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6 Appendix 

6.1 Eigenvalues for a Blossom 

Remark that an eigenvalue A and eigenvector v of matrix A satisfy Av = Xv, i.e. 

^V2 = Xvi (12) 
11, 

2^i-i + 2^^+i = ^ <i <n' + m (13) 

1 1 

"2^" 2^"+™"l ^ -^^n+m (14) 

Suppose A = cos(0) and v = (sin(0), sin(20), . . . , sin((?2 + m)0))"^. Then, by elementary 
trigonometric identities we note that ( fT2l) and ( fT3i) hold for every 0. On the other hand, (HM 
is equivalent to 

sin((n + m — 1)0) = sin(n0) + 2 cos(0) sin((?2 + m)0) 
which by using elementary trigonometric identities is equivalent to 

. ,'2n + m + l\ / m + 1 , 
sm (p cos (p = 0. 



2 V 2 

Therefore, is either 

2ki 2k2 - 1 

= o — ^ — °^ = — rr^ 

2n + m + 1 m + 1 

where ki and ^2 are arbitrary integers. Since cosine is a periodic function, it can be readily 
checked that cos(0i) attains all possible values over k = 1, 2, . . . , n + [m/2j and similarly for 
cos(02) over A; = 1, 2, ... , [m/2] . 
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6.2 Eigenvalues for a Bicycle 

If A is an eigenvalue of matrix A with eigenvector w, then we have 

1 



Xvi 
Xvi 

XVn^m+l 



2^2 - 2^^+i 



n + m + / — 1 



In the remainder, we separately consider two cases depending on whether either / or m 
is even, or otherwise. 
Case 1: / or m is odd. 

Without loss of generality, suppose / is odd. Let us introduce the following one-to-one 
linear transformation z = Sv where matrix S is defined by Zi = Vi+Vi-i+i, for z = 1, . . . , [//2J , 
and Zi = 2vi for i = [l/2\+l,...,n + l + m. It is not difficult to verify that S is non-singular 
and thus a matrix B such that A = S^^BS is similar to A and, therefore, has the same 
eigenvalues [8J [Theorem 1.3.3]. 

Using the transformation z = Sv and Av = Xv, we have 



and for i 



Xzi 

XZn+m+l 

[l/2\ and i = [l/2\ - 
Xzi 



-Z2 

2 

Z\l/2\ 
1 



1 



2 /t-r//t^t— J- 2 

1, ...,n + / + m — 1 
1 1 

-Zi-l + -Zi+l ■ 



Zn+l 



Notice that Xz = SAS 



-1; 



B£ and from the above identities 



B 



P 
Q R 



where P is a [//2] x [//2] tridiagonal matrix given by 



/ 

1/2 


V 



1/2 



1/2 




\ 


1/2 
/ 



Q is a (n + m + [//2j) X [//2] matrix with all elements equal to zero but the element in the 
first row and last column equal to 1/2, and R is a {n + m+ [l/2\) x (n + m+ [V2j) matrix 
that corresponds to a blossom with n+ 11/2] matched stem edges and m matched loop edges 
and is of the form ([9]) under Case 1 in Section 13.31 
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Using the properties of determinants of block matrices, we observe that eigenvalues of B 
consist of eigenvalues of matrices P and R. Therefore, the eigenvalues of matrix B, and by 
similarity of matrix A, are 

cos(^(|^), fc=l,...,r//2l, (15) 
cos 

COS ' 



A: = l,...,rm/21, (16) 
(2^)'^ = l'---'^+LV2j + LW2j (17) 



where (fT5|) are eigenvalues of matrix P, which is easily derived and thus omitted, and (fT6l) 
and f fT7|) are eigenvalues of R which we have already showed in Section 13.31 

It is not difficult to see that the above eigenvalues hold whenever either / or m is odd. 
Case 2: both I and m are even. 

We use a similar but different one-to-one transformation as under Case 1: Zi = + 
fori = 1, . . . ,1/2, Zi = Vi+Vi+i for i = 1/2+1, . . . , n+l+m/2, and Zi+n+i = Vn+i+i+Vm+n+i-i+i 
for 2 = 0, . . . , m/2. We have that 

>^Zi = ^Z2 

- _ 1 1 

— -Zl/2-l + 2^1/2 

_ 1 1 

XZn+l+m/2 — + 2^"+'+™/2+l 

, _ 1 

and for z = //2 + 1, . . . , n + Z + m/2 — 1 and i = n + I + m/2 + 1, . . . , n + m + Z — 1, 

, _ 1 1 

■^Zi — 2^i-i + 2^*+^ ■ 

Similarly as for Case 1, using the properties of determinants of block matrices, we have 
that the eigenvalues of A are 

cos A; = l,...,Z/2, 

cos — — , fc = 1, . . . , m/2, 

\ m + 1 / 

cos ( 7^ 777 ) , k = l,...,n + Z/2 + m/2- 1, 



n + 1/2 + m/2 
and — 1 

where (1, —1, 1, —1, . . . , 1, —1)^ is the eigenvector of eigenvalue —1. 
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